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INTRODUCTION 

This paper describes the multiphase effort in the development of tools 
for the analysis and design of two**element airfoil Systems — that is, air- 
foils with a slat or a flap at transonic speeds. The first phase involved 
the development of a method to compute the invlscld flow over such conflgura- 
tions. In the second phase the invlscld code was coupled to a boundary layer 
calculation program in order to compute the loss in performance due to 
viscous effects. In the third phase, which is currently in progress, an 
inverse code that constructs the airfoil system corresponding to a desired 
pressure distribution is to be developed. The symbols are defined in an 
appendix. 


TWO-ELEMENT AIRFOIL ANALYSIS 

The methodology used to compute the Invlscld flow over two-element airfoil 
systems has been described in previous reports (ref. 1-3). It will be 
summarized here only briefly. The flow field over such configurations is 
obtained by mapping the doubly connected physical domain into an annulus. 

The transformation Itself follows from the work of Ives (ref. 4), and it is 
the result of five sequential mappings. The first, a Karmar-Trefftz trans- 
formation, opens up the main airfoil into a near circle; the second turns the 
near circle into e perfect one. Then the second airfoil is opened up and 
transferred to the interior of the main circle. Finally, the second foil 
is mapped into a circle while managing to retain the shape of the first circle. 
With this transformation the entire physical space is mapped into the annular 
region between the two circles with infinity becoming a single point bet’ en 
them. 

The governing equations for the Inviscid, l^ . u tional compressible flow 
are written in this computational domain using the metric of the mapping. A 
potential function is then introduced to reduce the number of equations to one. 
The mapping introduces several singularities in the potential equation, but 
stace the mapping is analytic and the transformation is known everywhere, the 
sin^larltles can be removed analytically. The metric of the transformation 
is thus normalized with its value at the infinity point, and terms corresponding 
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to the free stream and the circulatory flows around each element are 
subtracted from the potential function. The result Is a second-order partial 
diff.<irentlal equation for a reduced potential function, whica is continuous 
and single-valued everywhere. As usual in transonic flow problems the equation 
is of mixed type and is solved numerically. The boundary conditions imposed 
are, of course, that the normal component of velocity be zero at each surface. 
The equation contains two lirculatlon constants that are evaluated by 
applying the Kutta condition at each trailing edge. 

In the computational domain a simple polar coordinate system emanating 
from the center of the annulus automatically generates an orthogonal grid 
in which both airfoils lie along grid lines. An accurate application of the 
boundary conditions is then made in a relatively straightforward manner. As 
described in reference 2, an additional stretching is used to concentrate 
points near the leading and trailing edges of both elements and to place both 
trailing edge points and the point of Infinity exactly at grid points. The 
numerical procedure employs standard relaxation techniques along with a non- 
conservative, type-dependent, rotated difference scheme. To make sure that 
the field is never swept at more than 90“ from the streamline direction, the 
computational domain is divided into four regions along the ring going through 
the infinity point, as shown in figure 1, and then the region over each 
of the airfoil surfaces is swept from the leading edge to the trailing edge. 
Figure 2 shows a typical computational grid as it appears in the physical 
plane. The appropriate high level of concentration of mesh points is obtained 
near the leading and trailing edges and in the gap. 

The configuration shown in figure 2 is a classic one, and the results of 
a computation for a Mach number of 0.6 and an angle of attack of 6* are given 
in Figure 3. Very large supersonic regions can be seen to be present on both 
elements. These are evident in figure 4, which gives the sonic line along 
with some computed streamlines. 

Transonic data for two-element airfoil systems are scarce, and there are 
little data at any speed where viscous effects are negligible. However, data 
recently made available by the David Taylor Research and Development Center 
have made possible the verification of the results of the method. The 
airfoil-slat combination shown in figure 5 was designed for lew-speed 
application on a circulation control wing. The unconventional back end of 
the main airfoil was designed to operate as a Coanda jet. A jet of high 
velocity air is ejected tangentially along the upper surface near the trailing 
edge of the main airfoil. The jet wraps around the rounded trailing edge 
entraining the outer flow. The result is an airfoil system with a very 
high circulation. The slat is deployed to prevent flow separation near the 
leading edge of the main conqjonent. The case shown in figure 5 is for a low 
Mach number, an angle of attack of 12“ , and zero blowing. Because of the 
small amount of aft loading, viscous effects on the main airfoil arc small 
and because of the high angle of attack there is very little separation on 
the slat. As a result there is very good agreement between the computations 
and the experimental data. Leading edge expansions are predicted correctly 
on both elements. The only discrepancy is on the lower surface of the slat 
where a small separation bubble is likely to exist. The lift coefficient on 
this configuration is 1.83. 
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DESIGN OF TVO-ELEMENT AIRFOIL SYSTEMS 


An analysis program can be of great help in the design of new airfoil 
configurations, but even so the analytical design could be a long process. 

The program would have to be run repeatedly with many contour modifications 
to achieve the desired flow characteristics. A computer program which could 
solve the inverse problem would indeed be valuable. An inverse program is 
one that generates the airfoil system that corresponds to a given pressure or 
velocity distribution. The development of such a program is currently under- 
way. The approach to the inverse problem is that described by Tranen in 
reference 6 for the single airfoil design. The method consists of prescribing 
a desired velocity distribution on a starting airfoil system. The input 
velocities are integrated on each of the two airfoil surfaces to get the 
potential on the boundaries of the computational domain. The Dlrlchlet 
problem in the annulus is then solved and the normal component of velocity 
on the boundaries is computed. The boundaries are not necessarily stream- 
lines now, but the computed normal velocities at the boundaries can be used 
to find nearby streamlines which are taken to be sought-after contours. This 
procedure is repeated until the normal velocity at the boundaries is driven 
to zero. But convergence is very fast. In the single airfoil design three 
or four iterations are usually sufficient. Of course, one could elect to 
design only one airfoil element or pieces of one or of both elements. In 
such cases a mixed boundary value problem needs to be solved, but this does 
not entail any additional numerical difficulties. 

A difficulty with airfoil design is that no full solution of the inverse 
problem has ever been given, even for the case of incompressible flow over a 
single airfoil. The problem arises because there is no assurance that a 
physically acceptable airfoil shape will follow from an arbitrarily prescribed 
pressure distribution. For a given surface pressure distribution an airfoil 
will generally not close at the trailing edge and the upper and lower surfaces 
may intersect. Provisions have to be made to allow the achievable pressure 
distribution to vary somewhat from the initially prescribed values to obtain 
closure. Also it is probably true that if we wish to design both airfoil 
elements completely one will not be able to specify the pressure distribution 
on one element without taking into account the effect on this element due to 
the flow characteristics required on the second element. But this situation 
shoul'’ not arise often in practice though. More likely will be the problem 
of redesigning only one of the airfoils or, even more commonly, only pieces 
of one element. This brings about the third type of design problem to be 
considered. In the cativ*} of a mixed design there are some regularity 
conditions that must be satisfied at the point where one switches from a 
geometry to a velocity boundary condition. This problem, as the trailing 
edge closure problem, is not unique to the two-element design, but is also 
present in the single airfoil design. It has not been dealt with satisfac- 
torily even in this case. 
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CONCLUDING REMARKS 


Beyond completion of the inverse code there are additional tasks that 
need to be undertaken to make the tools more useful to the designer* Melnik 
et al. have recently developed a more rational model for the viscous trailing 
edge* Its implementation in the analysis would certainly be an improvement* 
Also the viscous analysis code would have to be generalized to take into 
account the possibility of the boundary layers and/or wakes from the two 
elements merging as might be the case where the two elements are very closely 
spaced* In the design problem the desirable characteristics of a pressure 
distribution need to be determined* Clearly one wants shocks to be weak and 
to eliminate, or at least reduce, boundary layer separation* A bit harder to 
determine is what to design for, or try to achieve, on the surfaces along the 
slot between the two airfoils* Experimental data would serve as a useful 
guideline in this task* Shortcomings of previous designs and discrepancies 
between data and theory can give useful indications on which are the desirable 
pressure characteristics and which are the undesirable ones* 
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APPEIJDIX 


SYMBOLS 

A radius of point of infinity in annular donaln 
lift coefficient 
CP pressure coefficient 

blowing coefficient 
M^ free-stream Mach number 

r radial direction in annular domain 

R radius of inner ring corresponding to secondary airfoil in annular domain 
Re Reynolds number 

<x angle of attack 

6 azimuthal direction in annular domain 
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Figure A.- Computed streamlines for Clark Y airfoil with 30 percent 
Maxwell slat. M ■ 0.6: a ■ 6**. 
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Figure 5.- Computed and experimental surface pressures on a CCVJ airfoi; 
with 25“ slat. M - 0.1; a = 12“; C » 0. 
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Figure 7.“ Computed and experimental surface pressures on a CCW 
airfoil with 25® slat. M * 0.1; a = 4.6®; C = 0.1285. 
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Figure 8.- Computeci and experimental surface pressures on an 
NACA 64A406 airfoil with 7.8A slat. 
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